STABLE SELF-SIMILAR BLOW UP FOR ENERGY SUBCRITICAL WAVE 

EQUATIONS 



ROLAND DONNINGER AND BIRGIT SCHORKHUBER 



Abstract. We consider the semilinear wave equation 

dtip - Aip = i^r _i v 

for 1 < p < 3 with radial data in R 3 . This equation admits an explicit spatially homogeneous blow 
up solution if) T given by 

where T > and k p is a p-dependent constant. We prove that the blow up described by ip T is 
stable against small perturbations in the energy topology. This complements previous results by 
Merle and Zaag. 

1. Introduction 

We study the Cauchy problem 

/ d<v - a?/, = ivr v n n 

\ #)] = (/,</) { } 

for tp : [0, oo) x R 3 — > M. and 1 < p < 3. The associated energy reads 

E(^(t,-),Mtr)) = \ [ {\Mt,x)\ 2 + \V x ij(t,x)\ 2 )dx--^— [ m,x)\ p+1 dx. 

2 JR3 p + 1 J R 3 

Furthermore, Eq. (jl.ip is invariant under the scaling transformation 

2 

ip(t, x) ^ ip x (t, x) := A p- 1 ip(t/\,x/\) 
for A > and the energy scales according to 

EtyxQ, •), d t Mt, •)) = \-*=*E(i/,{t/\, -),A(t/X, •))• 

Thus, Eq. (jl.ip is energy subcritical if 1 < p < 5, critical if p = 5 and supercritical if p > 5. 
Note, however, that the usual blow up heuristics (energy conservation prevents the solution from 
shrinking to ever smaller scales in the subcritical case) do not apply here since the energy is not 
positive definite. 

Local well-posedness of Eq. (11. ip in ^T 1 (]R 3 ) x L 2 (R 3 ) follows by standard arguments (use 
Duhamel's formula and Sobolev embedding to set up a contraction scheme). However, Eq. (jl.lj) is 
not globally well-posed and it is well-known that initial data with negative energy lead to singu- 
larity formation in finite time [10]. More explicitly, one may look for self-similar solutions which 
are by definition invariant under the natural scaling. A particularly simple example of this type is 
obtained by neglecting the Laplacian altogether and solving the remaining ODE in t which yields 
the solution 

iP T (t,x) = Kf T (T-t)-^ 
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where T > and kq = |~z7[p • We refer to as the fundamental self-similar solution. Although ijj T 
is homogeneous in space, by finite speed of propagation it can still be used to construct compactly 
supported smooth initial data that lead to blow up at time t = T. In fact, one expects that Eq. 
admits many self-similar blow up solutions, even in the radial case. At least for p = 3 this was 
proved by Bizoh et. al. pQ (there are similar results for p > 7 |3j). We remark that the situation 
for the corresponding problem in one space dimension is fundamentally different since in this case 
there exists a unique (up to symmetries) self-similar blow up solution. 

In order to understand the dynamics of Eq. (jl.ip it is important to analyse possible blow up 
scenarios. In two remarkable papers |12| Merle and Zaag proved that any blow up solution 
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of Eq. (jl.ip blows up at the rate (T — t) p- 1 . However, the precise shape of the blow up profile 
depends on the data. Numerical work by Bizoh, Chmaj and Tabor [2] suggests that the blow up 
described by ip T is the "generic" one. To be more precise, they observe in their simulations that 
the future development of "generically" chosen radial blow up data converges to ifj T locally near 
the center of spherical symmetry. In the present paper we analyse the stability of and obtain 
the following result, see Theorem 11.21 below for the precise formulation. 

Theorem 1.1 (main result, qualitative formulation). There exists an open set (in the energy 
topology) of radial data that lead to blow up via ip T . In other words, the blow up described by tp T is 
stable. 

Here, "energy topology" refers to the topology generated by the energy of the free wave equation. 
In fact, we use a slight modification thereof which has better local behavior, see below. Before 
going into more details, let us briefly comment on possible extensions of our result. First of all, 
the restriction to radial data is only technical and can be quite easily removed. However, since the 
nonradial case is not expected to reveal any new interesting phenomena, we stick to radial data 
so as not to obscure the exposition by too many technicalities. A more important aspect concerns 
the admissible values of p. It should be possible to extend our result to the full subcritical regime 
1 < p < 5 with some modifications. For instance, it is necessary to require slightly more regularity 
(e.g. of Strichartz type) if 3 < p < 5 in order to control the nonlinearity. This is reminiscent of 
the standard well-posedness theory for semilinear wave equations and we will pursue this matter 
elsewhere. We even claim that the stability of ip T holds true in the full supercritical regime p > 5 
albeit in a topology significantly stronger than the energy, cf. [7J, [6] for analogous results on 
supercritical wave maps. 

In a certain sense our work is complementary to the results of Merle and Zaag [12J. In their 
approach they consider any blow up solution whereas our method is perturbative near i(j t . Conse- 
quently, the strength of their result lies in its generality: it is a statement about all possible blow 
up scenarios. However, it is not true that every blow up solution converges to ip T since there exist 
many self-similar profiles. We remark that the situation in the one dimensional case is entirely dif- 
ferent due to the uniqueness of the self-similar solution. As a consequence, the Merle-Zaag method 
even yields profile convergence for the problem in R 1+1 |13]. Furthermore, a number of beautiful 
results emerged from this approach, see, e.g. [15], [H], [16], [4]. 

Our method of proof follows essentially our work on the wave maps problem [7J, [6]. However, 
there are important differences. First of all, the problem at hand is energy subcritical and this 
allows us to work in the energy topology whereas for the wave maps problem we had to require 
more regularity. Furthermore, the present result is completely rigorous and there is no need for 
any numerical input as was the case in |6j. This is because the nonself- adjoint spectral problem 
associated to the linear stability of ip T can be solved explicitly in terms of hypergeometric functions. 
Finally, we have improved the method in order to avoid one additional degree of differentiability 
which was still necessary in [6j. Consequently, the motivation for this paper is in fact twofold: 
First, our result is intended to complement the work of Merle and Zaag in order to obtain a fairly 
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satisfactory description of the blow up behavior of Eq. (jl.ip . Second, we wanted to demonstrate 
the wide applicability of our methods developed in [7] and [6]. 

Our approach is functional analytic. We first introduce new coordinates adapted to self-similarity 
and linearize the equation around the solution tp T . The resulting linear problem involves a highly 
nonself- adjoint operator (due to the change of coordinates) and it is therefore studied by semigroup 
methods. This yields the linear stability of if 1 '. In fact, this stability is modulo the time translation 
symmetry of the problem which manifests itself in the form of an unstable mode of the linearized 
operator. We remove this instability by a Riesz projection. The nonlinear stability is then proved by 
a fixed point argument where we have to take into account the instability caused by the symmetry. 
In order to deal with this problem we employ an infinite-dimensional version of the Lyapunov- 
Perron method from dynamical systems theory: we force nonlinear stability by modifying the data. 
In a last step we undo this modification by shifting the blow up time. 

1.1. Formulation of the Cauchy problem and the main result. We restrict ourselves to 
radial solutions and study the initial value problem in the backward lightcone Ct of the blow up 
point (T, 0) which is defined by 

C T :={{t,r):te(Q,T),re[0,T-t]}. 

More precisely, we consider 

' Mt,r) -Ar{t,r)-^Mt,r)-\iP{t,r)\P~ l ^(t,r) = for (t,r)eC T n „s 

ip(Q,r) = f(r),ik(0,r) = g(r) for r € [0, T] 1 ' > 

with given initial data (f,g)- We are interested in the stability of ip T under small perturbations tp. 
Thus, we insert the ansatz ip = ip T + ip into Eq. (|1.2j) and expand the nonlinearity according to 

where Nt denotes the nonlinear remainder. Since ip T (t,r) > for all (t,r) <E C T we obtain 



Ptt ~ Prr ~ ~<Pr ~ P{^ T ) P V ~ N T (<f) = in C T 

<p(P,r) =/(r)-V T (0,r),(^ f (0,r) =g(r)-^(0,r) for r G [0, T] 



1.2. Energy norm. We want to study the Cauchy problem in a backward lightcone and, since our 
approach is perturbative, we need a local energy norm derived from the conserved energy associated 
to the free equation 

2 

Pit ~ Vrr <Pr = 0, (1.4) 

r 

which is given by 

f oo 

r 2 [ip t {t,rf + ip r (t,r) 2 ]dr. 

o 

However, this does not define such a norm due to the lack of a boundary condition for ip at r = 0. 
By integration by parts (and assuming sufficient decay at infinity) it can be easily seen that the 
above expression for the energy is equivalent to 

POO 

E(<p)= r 2 ip t {t,rf + [rp r (t,r) + p(t,r)] 2 dr. (1.5) 
Jo 

Another way of motivating this is to define (p := rep, such that Eq. (jl.4|) transforms to the 1 + 1 
wave equation 

<Ptt - ftrr = 

with conserved energy 

Pt(t,r) 2 + <p r (t,r) 2 dr. 



Writing this expression in terms of the original field yields (|1.5p . For (/,<?) G C 1 ^,/?] x C[0,R], 
R > we define 

Kfi9)\\e(R) '•= J \rf\r)+f(r)\ 2 dr + ^ r 2 \g(r)\ 2 dr. (1.6) 

r f'( r ) + = implies f(r) = p and the requirement / G C [0, i2] yields c = such that || • \\gfm 
defines a norm on C 1 [0, R] x C[0, i?]. We consider the expression (|1.6|) in the backward lightcone 
of the blow up point (T, 0) and insert the fundamental self-similar solution to obtain 

\\^ T (t,-),^(t,-))\\ s{T _ t) =C p {T-ty&% (1.7) 

where C p > denotes a p-dependent constant. Evidently, as t — > T— , this quantity blows up in 
the energy subcritical case, i.e., for 1 < p < 5. 

Theorem 1.2 (main result, quantitative version). Fix 1 < p < 3 and e > 0. Let (f,g) be radial 
initial data with 

||(/, 5 )-(^(o,-),^(o,-))IU(|) 

sufficiently small. Then there exists a T > c/ose to 1 such that the Cauchy problem 

#)] = (/,</) 1 j 

/ias a unique radial solution tp : Ct — > K which satisfies 

{ t - t)^\\m, -),Mt, •)) - (^ T (t, (*, -))iU(T-t) < c £ (r - t)W-« 

/or a// i G [0, T) where uj p := max |— 1, | — ^rx| and C £ > is a constant which depends on e. 

2. Transformation to a first order system and similarity coordinates 

2.1. First— order formulation. By setting <p(t,r) := rip(t,r), Eq. (|1.3p transforms into 

<Ptt ~ ftrr ~ p(ip T ) p ~ 1 <p - rN T (£) = in C T (2 ^ 

<^(0,r)=r[/(r)-V T (0,r)],^(0,r)=r[ 5 (r)-V t T (0,r)] for r G [0,T] 1 ' J 

together with the boundary condition <p(t, 0) = for all t. We rewrite Eq. (12, ip as a first-order 
system by introducing new variables 

2 2 

</?l = (T-t)p- 1 ^, 992 = (T- t)^- 1 ^. 
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Thanks to the boundary condition we have <p(t,r) = (T — t) p- 1 L (f2(t, r')dr' . Note that the 
nonlinearity transforms according to 

N T ((T - ty^r" 1 / p 2 ) = (T — ty^Nir' 1 J y? 2 ) 

where J y>% is shorthand for J Q r (f2(t, r')dr' and 

N(x) = + xf 1 ^^ 1 +x)- k^ 1 - P k x. (2.2) 

Thus, Eq. (I23D reads 

3t¥>i =9 r ^ 2 -^(T-t)-Vi+^o(T-t)- 2 /V2 + r(T-t)- 2 iV(r- 1 / ^ 2 ) 



d t p 2 = d r <Pi - ^j(T - t) V2 

^(0,0 =T^r [ 5 (r)-^ T (0,r)] 
[ p 2 (0,r) =T^T[r/ / (r) + /(r)-V T (0,r)] , 

1 



in Ct 



for r G [0, T] 



2.2. Similarity coordinates. We transform the system to similarity coordinates (r, p), which are 
given by 

r = -io g (T-t), p = j^- t - 

The cone Ct gets mapped to the infinite cylinder := {(r, p) : r > —logT,p £ [0,1]} and by 
setting 

4>j{i~,p) := (fj(T - e~ T ,e' T p) 

for j = 1 , 2 we obtain 



d T 4>i = -pd p (pi + 9p0 2 - ^rr^i + P«o J P 02 (t, s)ds + piV(p 1 f Q p 2 (r, s)ds) 



9 T 02 = -pd p 4>2 + 5 p 0i - ^j- 

r D+i 

0i(- log T,p) = p 



in Zt 



p+i 

T*-ig{Tp) 



i 

p-i 



(2.3) 



i 

.p-i 



for p £ [0, 1]. 



M~ log T, p) = T—i [Tpf'(Tp) + /(Tp)] - «*" 
Tracing back the above transformations, the original field tp as well as its time-derivative can be 



reconstructed according to 



ip{t, r) = tp T (t, r) + (T — t) p-ir 



-i 



6 2 (-log(T-t),^)dr', 



^(t, r) = $(t, r) + {T- ty^r^M- log(T - i), ^). 



(2.4) 



Note that most of the expressions we are going to define below will depend on p. However, for 
the sake of readability we will not indicate this dependence explicitly, but consider p to be fixed, 
where we restrict ourselves to 1 < p < 3. 

3. Linear Perturbation Theory 
In this section we consider the linearized problem 

dr4>i = -pd p (f)i + d p 4>2 - ^=4 01 + PK J p (f>2 (t, s)ds 



d T 4>2 = -pd p 4>2 + dptpi - ^4>2 



l(- log T,p) = p 



p+i 



h(~ log T,p) = T 2 /(p-D [Tpf'(Tp) + f(Tp)] 




(3.1) 



which has already been studied in [5]. Nevertheless, in order to present a consistent picture we 
summarize known results and supplement them by some new aspects which will be important for 
the nonlinear theory (see Lemma |3.7[) . 

3.1. Well-posedness of the linearized equation. Let H := L 2 (0, 1) x L 2 (0, 1) with the usual 
inner product. We define operators (L ,V(L )) and V G B(J-L) by 

V[L Q ) := {u G C^O, 1] x C^O, 1] : m(0) = 0}, 



L u(p) : = 



u' 2 (p) ~ pu[(p) - ^lMp) 
u'i(p) - pu' 2 {p) ~ ^iMp) 



and 



L'u(p) := 



pK f£ii2(s)ds 



where u = {ui,U2) T . It is easy to see that V is a compact operator, which will play an important 
role later on. 



Lemma 3.1. The operator Lq is closable and its closure Lq generates a strongly continuous one- 
parameter semigroup Sq : [0,oo) — > B{T~L) satisfying ||So(t)|| < ef ,pT for all r > and oj p := \ — -^y. 

Proof. The claim is a consequence of the Lumer-Phillips Theorem (see [8], p. 83, Theorem 3.15). 
Indeed, a simple integration by parts yields the estimate 

Re(Z u|u) < (| - ^) ||u|| 2 

and 7=j — — j < 0. Furthermore, for A := 1 — > uj p the range of A — Lq is dense in %. This follows 

from the very same calculation as in the proof of Lemma 2 in [5]. Since Lq is densely defined, the 
Lumer-Phillips Theorem applies. □ 

Corollary 3.2. The spectrum of Lq is contained in a left half plane, 

o-(Lq) C {A G C : ReX < u p } , 

with ojp = \ — rzj and the resolvent of Lq satisfies 

for all A G C with ReX > Cj p . 

Proof. The structure of the spectrum as well as the resolvent estimate follow by standard results 
of semigroup theory (see [8], p. 55, Theorem 1.10). □ 

The next corollary is a consequence of the Bounded Perturbation Theorem (see [8], p. 158). 

Corollary 3.3 (well-posedness of the linearized equation). The operator L := Lq + L' , T>(L) := 
T>(Lq) generates a strongly continuous one-parameter semigroup S : [0,oo) —> B(T-L) satisfying 

\\S(t)\\ < e (^+P K 0)r 

for all t > and Cj p = \ In particular, the Cauchy problem 

£*(r) =L$(r) for r > — log T 
$(-logT) = u 

/ias a unique solution given by 

$(r) = 5(r + logT)u 

/or u G T>(Lq) and all r > — log T. 

3.2. Properties of the generator. We obtain a more explicit characterization of T>(L) in order 
to be able to describe the spectrum of L. 

Lemma 3.4. Let u G V(L). Then u G C[0, 1) x C[0, 1) and ui(0) = 0. Furthermore, for f e H 
the equation (A — L)u = f implies 



rp rp 

ui(p) = pu 2 (p) + (a + - 1) y u 2 (s)ds - y^ 



f 2 (s)ds 



and 

-(1 - pV(p) + 2 (A + ^) pu'(p) + ((A + (A + ^ - l) - p Ko ) u(p) 

p (3 2) 

= + P/ 2 (P) + (A + ^i) J" f 2 (s)ds 

in a weak sense, where u G H? oc (0, 1) n C[0, 1] n C^O, 1) is denned oy u(p) := j^u 2 (s)ds. 
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Proof. Let u G T>{L) = V{Lq). By definition there exists a sequence (uj) C V{Lq) C C^fO, 1] x 
C 1 [0, 1] such that \ij — > u and Lo u j — > ^o u i n By combining the expressions for the individual 
components in an appropriate way we infer that (1 — p 2 )u'ij and (1 — p 2 )u' 2 j are convergent sequences 
in L 2 (0, 1). Thus u\,u 2 G H 1 ^, 1— e) >■ C[0, 1— e] for any e G (0, 1). This guarantees the boundary 
condition «i(0) = 0. 

Let f G T~L and A G C. Then (A — L)u = f implies 

p 

(A + ^i)ui(p) +pu' 1 (p) - u' 2 (p) - pK f u 2 (s)ds = h(p) 

v 

(A + ^i)Mp) + Pu' 2 (p) ~ u'i{p) = f2(p) 

in a weak sense. Thanks to the boundary condition we obtain from the second equation that 

rp rp 
ui(p) = pu 2 {p) + (A + ^rr - 1) J u 2 (s)ds-J f 2 (s)ds. 

Inserting this into the first equation yields 

-(1 - p 2 )u' 2 (p) + 2 (A + ^) pu 2 {p) + ((A + ^) (A + ^ - l) - pK ) /n 2 ( S )d S 

= flip) + p/ 2 (p) + (A + ^t) / h{s)ds. 

We set := f£ u 2 (s)ds and obtain Eq. ([372]) . Finally, u 2 G L 2 (0, 1) implies u G ^(O, 1) 
C[0, 1] and u G V{L) yields u G H 2 QC {0, 1) D C^O, 1). □ 

In order to improve the rough growth estimate given in Corollary 13.31 we analyse the spectrum 
of the generator. The next two Lemmas characterize the spectral properties of the generator L 
sufficiently accurate. 

Lemma 3.5. For the spectrum o~(L) of the generator L we have 

o-(L) C {A G C : ReX < max{w p , -1}} U {1} 

where uj p = ^ — . 

Proof. Set M := {A G C : ReX < max{w p , -1}} U {1}. Let A G a(L). If ReX < u p then A G M. 
So let us assume that ReX > uj p . Then, by Corollary 13.21 A G cr(L) \ o-(Lq) and the identity 
A — L = [1 — L' Ri (X)](X — L ) together with the spectral theorem for compact operators imply 
that A G o~p(L). Thus, there exists a nontrivial u G F){L) such that (A — L)u = 0. By Lemma 
13.41 this implies the existence of a weak solution u of Eq. (|3.2|) with right hand side equal to zero. 
Recall that u G i/ 1 (0, 1) and u(0) = 0. We transform Eq. (13. 2p by substituting p t- > z := p 2 to 
obtain the hyper geometric equation (recall that kq = z^zjp ) 

z(l - z)v"(z) + [c - (a + b + l)z]v'{z) - abv(z) = (3.3) 

where v (z) := u{^/z) and the parameters are given by a = ^(A — 2), b = ^(A + jpf), c=\- For A ^ 
1 — a fundamental system around z = 1 is given by = 2-^1(0, 6; a + 6 + 1 — c; 1 — 

and i)i(z) = (1 — z) c ~ a ~ h 2 F\(c— a, c— b; c+ 1 — a — b; 1 — z) where 2 F\ is the standard hypergeometric 
function, see e.g., [IT]. The exponent c — a — 6 = 1 — -^j — A vanishes for A = 1 — ^-j- and in this 
case one solution is still given by v\ and the second one diverges logarithmically for z — )■ 1. Since 
we assume ReX > \ — ^y, v must be a multiple of v\ for u to be in i? 1 (0, 1). Around z = there 
is a fundamental system given by {vo,vq}, where vq(z) := z l ~ c 2 F\{a + 1 — c, b + 1 — c; 2 — c; z) 
and i)o(<2) := 2^1 («, c; z). Thus, there exist constants ci, c 2 such that = c\Vq + C2t>o- In order 
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to satisfy the boundary condition u(0) = 0, the coefficient c\, which can be given in terms of the 
Gamma function [17J 

F(a + b + 1 - c)r(l - c) 



Cl r(a + l-c)r(6 + l-c)' 
must vanish. Consequently, c\ = if and only if a + 1 — c or 6+1 — c is a pole, which yields 
^ (A — 1) = —k or ^ + 2±i = — k for a k G No- This implies that A is real and 

AG{wel:w>w p A(w = l-2fcVw = -2k - , fc = 0, 1, . . . } . 

Since 1 < p < 3 we have Cj p < — ^. If w p > — 1 then the only possibility is A = 1. If Cj p < — 1 then 
either A = lorA<— 1. In any case we conclude that A G M. □ 

Lemma 3.6. The eigenvalue 1 G o~p(L) has geometric multiplicity equal to one. The associated 
geometric eigenspace is spanned by 

g(p) : = ( S P ) • ( 3 - 4 ) 
In £/ie following g 6e referred to as the symmetry mode. 

Proof. Note that g G D(L) and a straightforward calculation yields (1 — L)g = 0. In particular by 
Lemma |3, 41 and the definition of kq we infer that 

CP 

9i(p) = P92(p) + 92(s)ds (3.5) 

and 

- (i - p 2 W'(p) + 2J ^?p9'(p) - 2 -^g{ P ) = o (3.6) 

for g(p) := f£ g2(s)ds = p. Suppose there is another eigenfunction g for A = 1. Then this 
corresponds to another (weak) solution g(p) := f£ g~2(s)ds of Eq. (|3.6p . A fundamental system of 
Eq. (|3.6p is given by {ho, hi}, where ho(p) = p and 

h 1 {p) = {l-p 2 )-^h l (p) 

with/ii(p) = 2^1(1, l-jpr; \]P 2 ) and^i(l) ^ for 1< p < 3. However, by Lemma E31 5 G C[0,1] 

and thus it must be a multiple of ho = g. Therefore, there exists a constant c G C such that 

p rp 
g 2 (s)ds = c g 2 {s)ds 
J 

and we infer that g 2 = cg 2 - Eq. (|3.5p implies 51 = cg\ and we conclude that g = eg. □ 

3.3. Spectral projection and linear time evolution restricted to the stable subspace. 

The symmetry mode is an explicit example of an exponentially growing solution of the linearized 
equation. However, its origin will only become clear in the course of the nonlinear perturbation 
theory. We will see that it is due to the time translation invariance of the problem and thus, we do 
not consider this instability as a "physical" one. The aim of this section is to remove the symmetry 
eigenvalue A = 1 via a Riesz projection and to obtain a growth estimate for the solution of the 
linearized equation on the stable subspace. We define a projection operator 

P = ^ J R L (X)dX, (3.7) 

where 7 is a circle that lies entirely in p(L) and encloses the eigenvalue 1 in such a way that no other 
spectral points of L lie inside 7. The projection P commutes with L in the sense that PL C LP 
and as a consequence, P also commutes with the semigroup generated by L, i.e., PS(t) = S{t)P 
for r > 0. We define subspaces M. = Prl and N = (1 — P)TL which decompose the operator L into 



parts living on Ai and M, respectively. Let Lj\f be denned by Lj^u := Lu with V(Lj^) = V{L) C\N 
(Lm is then defined analogously) . Since M and Ai are closed subspaces we can regard Lj\f and L_A4 
as linear operators on the Hilbert spaces Af and Ai, respectively, with spectra (j(L_m) = {1} and 
a(Lj\r) = a(L)\{l}. In the following we call Ai the "unstable subspace". The operator L is not 
self-adjoint and therefore, the next result is nontrivial and crucial for the nonlinear perturbation 
theory. 

Lemma 3.7. The unstable subspace Ai is spanned by the symmetry mode, i.e., P% = (g) and the 
algebraic multiplicity of 1 £ o~ p (L) is one. 

Proof. The case dim Ai = oo can be ruled out by an abstract argument: if dim.A/f = oo then, by [9J, 
p. 239, Theorem 5.28, 1 would belong to the essential spectrum of L which is stable under compact 
perturbations (see [9J p. 244, Theorem 5.35). However, 1 o-(Lq) and this yields a contradiction. 
We conclude that Lm is in fact a finite-dimensional operator. Since 1 is an eigenvalue of and, 
according to Lemma [3. 6 1 the corresponding geometric eigenspace is spanned by g, we obtain g £ Ai 
and thus, (g) CM. 

It remains to prove the reverse inclusion. Note that (1 — Lj^) is nilpotent since is the only 
eigenvalue, i.e., there exists an m £ N such that (1 — Lj^) m \i = for arbitrary u £ Ai. If m = 1 
then At C ker(l — Ljm) = (g) and we are done. Suppose that m > 2. Then there exists a nontrivial 
v £ rg(l — Ljm) such that (1 — Lj^)v = 0, i.e., v £ ker(l — Lj^) and v must therefore be a multiple 
of the symmetry mode. This shows that there exists a u £ T>{L_m) with (1 — Ljk)u = eg. We will 
show that this leads to a contradiction. We set c = 1 without loss of generality. Suppose there 
exists a function u in T>(L) such that (1 — L)u = g. Then, by Lemma 13.41 

- (1 - pV(p) + 2J ^? P AP) - 2J ^<P) = 9(P) (3.8) 

for u(p) := f£u 2 (s)ds and g(p) := gi(p) + pg 2 (p) + |rr /o 92(s)ds = ^-p. For the homogeneous 
equation we have the fundamental system {ho, hi} introduced in the proof of Lemma 13.61 where 

2 ~ ~ ~ 

ho(p) = p and h\(p) = (1 — p 2 ) p^hiip) with hi continuous on [0, 1] and hi(0) / 0. The Wronskian 
is given by 

w{h ,hi) = -{\- p 2 )"^ 1 

and thus, a solution of the inhomogeneous equation must be of the form 

fP 2 _j_ fP 2 2 

u(p) = c h (p) + cihi(p) - h (p) / hi(s)g(s)(l - s )p~ 1 ds + hi(p) / h {s)g(s)(l - s ) p- 1 ds 

J po J pi 

for some constants cq,ci £ C and po,pi £ [0,1]. The boundary condition at p = implies 
ci = — r° h(j(s)g(s)(l — s 2 )~ds and inserting the definitions of /iq, hi and g yields 



'Pi 

H> . „ , , 2 - /'/' 



u(p) = c 0P - 3 -S+±pJ shi(s)ds + ^(l- p 2 )- — hi(p) J o s 2 (l-s 2 ) — ds. 

Since u belongs to C[0, 1] (Lemma 13. 4p . we must have J 1 s 2 (l — s 2 )p~ T ds = 0. However, this is 
impossible since the integrand is strictly positive for all s £ (0, 1). □ 

In order to improve the growth estimate in Lemma 13.31 we apply a well-known theorem by 
Gearhart, Priiss and Greiner. To this end we need the following result, which states that the 
resolvent is uniformly bounded in some right half plane. In the following we set H a := {A £ C : 
ReA > a} for a £ R. 

Lemma 3.8. For any e > there exist constants ci,oi > such that 

\\Rl(X)\\ <ci 

9 



for all A G Ha, p+£ with |A| > c 2 . 

Proof. Fix e > and let A G Hu p +e where A {1 — 1}. We use the identity 

#l(A) = i? Lo (A)[l-L / J R Lo (A)]- 1 
to obtain uniform bounds on the resolvent for |A| large. By definition of U we have 

L>Raw=( pKoK[R Q° m2 ) 

where K : L 2 (0, 1) — > L 2 (0, 1) is defined by Ku(p) = J* u(s)ds. For f E 7i consider the equation 
(A — Lq)u = f. Its solution is given by u = Rl (\){. Lemma [33] yields 

[i? i0 (A)f] x (p) = (A - 1 + ^)K[i? Lo (A)f] 2 (p) + p[R Lo (\)fUp) - Kf 2 (p). 

The estimate in Lemma 13.21 implies 

llfll 



[^o(A)f],-||x 2 (o )13 < ll^ (A)f|| < - |ReA _^, 



for j = 1,2 and we obtain 



|^[ J R Lo (A)f] 2 || L2(0jl) < 



IA-i + ^tI' 

Thus, for |A| sufficiently large, the Neumann series 

oo 

[l-L'R I ^(X)]- 1 = y £ i [L'R Lo (\)] k 

converges and the claim follows. □ 

We conclude the linear perturbation theory with an estimate of the linear evolution on the stable 
subspace. 



Proposition 3.9. Let P be the spectral projection defined in Eq. \3. 7| ) and set 

u p :=max{-l,i- 

Then, for any e > 0, there exists a constant C £ > such that the semigroup S(r) given in Corollary 
HO satisfies 

\\S(t)(1 - P)f\\ < C £ e^\ +e ^\\{\ - P)f\\ (3.9) 
for all t > and f € H. Furthermore, S{r)Pf = e T Pf. 

Proof. The operator Lj\f is the generator of the subspace semigroup S(t)\j^ and its resolvent is 
given by Rl(X)\j^. The first estimate follows from the uniform boundedness of the resolvent in 
-ff^+e (Lemma 13 8j) and the theorem by Gearhart, Priiss and Greiner (see for example [8J, p. 302, 
Theorem 1.11). The second assertion follows from PH. = (g) and the fact that g is an eigenfunction 
of the linear operator L with eigenvalue 1. □ 
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4. Nonlinear perturbation theory 

4.1. Preliminaries. Now we turn to the full nonlinear problem. The following two lemmas will 
be used frequently. 

Lemma 4.1. If u G L 2 (0, 1) i/ien u, defined by u(p) := f£ u(s)ds, belongs to L°°(0, 1) and 
satisfies 

IHIl°°(o,i) ^ IMIl 2 (o,i)- 

Proof. First note that p i— > f£ u(s)ds is a continuous function on [0, 1] for it G L 2 (0, 1). Using the 
Cauchy-Schwarz inequality we esimate 



Hp)\ 



i p 

— fu{s)ds 



< u 



L 2 (0,1) 



for p G (0, 1]. Taking the essential supremum yields the claim. □ 

We will also use Hardy's inequality in the following form. 
Lemma 4.2. For u G L 2 (0, 1) we have 

- 1 | f£u(s)ds 



P 2 



4.2. Estimates for the nonlinearity. From now on we restrict ourselves to real-valued functions. 
We introduce a function n : R x [0, 1] — > R defined by 

/ ^_ ^_ _e_ 

n(x, p) := pi + + x) — pK x - k^ 1 

cf. Eq. (|2.2p . It is easy to see that 

i ^ m <r / /°l x l 2 \ x \ < 1 

K^p)I<| Hx |p W > L 

A convenient way to write this is \n(x,p)\ < p\x\ 2 (x) p ~ 2 with the "japanese bracket" (x) := 
\/l + |x| 2 . In the following we denote by B\ and B\ the open unit balls in L 2 (0, 1) and H, re- 
spectively. To (formally) define the nonlinearity we introduce an operator A : L 2 (0, 1) — > L 2 (0, 1), 

1 [ p 

Au(p) := - / u(s)ds. 
P Jo 

An application of Hardy's inequality shows that A is bounded. We set 

N(u)(p) :=n(Au(p),p). 

Lemma 4.3. The operator N maps L 2 (0, 1) into L 2 (0, 1). Furthermore, there exist constants 
c\ , C2 > such that for u, v G B\ 

\\N(u)\\ L2 < Cl \\uf L2 

and 

\\N(u) - N{v)\\ L 2 < c 2 (\\u\\ L 2 + \\v\\ L2 )\\u - v\\ L 2. 

Proof. Note that for 1 < p < 3 the function n defined as above is at least once continuously 
differentiable with respect to x and we have the bound 



\din(x,p)\ < p\x\(x) 
li 



p-2 



for all x G R and /) G [0, 1] which, in particular, implies t?in(0, p) = and hence, N(0) = 0. By the 
fundamental theorem of calculus we infer that for x, y G R 

Kx,p) - n(y,p)| < \x-y\ \din(y + h(x - y) , p)\dh 

Jo 

^ P\ x ~y\ \y + h(x - y)\(y + h(x - y)) p ~ 2 dh 



< 



i _ i J M + M pe (1,2] 
p|x y| \ |x|<x>f- 2 + |yKy>P- 2 p 6(2,3] 



Now we prove the estimate for the nonlinear operator N. The following argument works only 
for 1 < p < 3, since for higher exponents the singular factors at p = can no longer be controlled. 
For u, v G L 2 (0, 1) we write u(p) := f£ u(s)ds and S(p) := f^vfyds. We distinguish two cases. If 
p G (1,2] we readily estimate 



\\N(u)-N(v)\\ 2 L2 = [ 1 \n(Au(p),p)-n(Av(p),p)\ 2 dp 



o 



< 



1 pVu(p) - Av{p)\ 2 {\Au{p)\ 2 + \Av{p)\ 2 )dp 



o 



£ (Nl£« + / |^(p) + ^(/5)| 2 rfp 

JO 

^(Nlis + Nli*-* 

by Lemma 14.11 and Hardy's inequality. On the other hand, if p G (2,3], we proceed similarly and 
obtain 

\\N(u) - N(v)\\ 2 L2 < f 1 p 2 \Au(p) - Av(p)\ 2 (\Au(p)\ 2 (Au(p)) 2 ^ + \Av(p)\ 2 (Av(p)) 2 ^)dp 
Jo 

< I'' p 3 -P\Au(p) - Av(p)\ 2 



x (\ P -\u{p)\ 2 {p-^u{p)) 2{ v-V + \p-k(p)\ 2 (p-k( P )) 2( p- 2 ^ dp 

< ( MU\\u\\^) 2{p - 2) + IMlLdMb) 2 ^) ll- - v\\ 2 L2 



again by Lemma l4.1l and Hardy's inequality. Since N(0) = we immediately conclude the bound- 
edness of N on L 2 (0, 1). In particular, we have ||iV(u)||x,2 < IMlia for u G B\. For u,v G B\ the 
above estimates yield 

\\N{u) - N{v)\\ L 2 < (\\u\\ L 2 + \\v\\ L 2)\\u - v\\ L 2 
as claimed. □ 



Finally for u = (ui,U2) T G % we define the vector valued nonlinearity by 

N(u) : = 



N{u 2 ) 


Lemma 4.4. The nonlinearity N maps T~L into %, N(0) = and there exist constants c±,C2 > 
such that for u,v G B\ 

||N(u)|| < ci||u|| 2 

and 

||N(u)-N(v)|| <e 2 (||u|| + ||v||)||u-v||. 
Furthermore, N is Frechet differentiable at and DN(0) = 0. 
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Proof. For u,v£Bi we apply the result of Lemma 14.31 to obtain 

||N(u) - N(v)|| 2 = \\N(u 2 ) - N(v 2 )f L2 < (\\u 2 \\l 2 + |hll| 2 )N - V2\\h < (IN 3 + l|v|| 2 )||u - v|| 2 . 
This implies 

||N(u) - N(v)|| < (||u|| 2 + ||v|| 2 )l||u - v|| < (||u|| + ||v||)||u - v||. 

We have N(0) = which implies ||N(v)|| < ||v|| 2 . In particular, there exists a constant c indepen- 
dent of v such that 

l|N(v)|| 

— n — n — — c v • 
||v|| 

Since the left hand side vanishes in the limit v — > 0, we infer that N is Frechet differentiable at 
zero with DN(0)=0. □ 

4.3. Abstract formulation of the nonlinear equation. We turn to the full nonlinear problem 
and write Eq. (|2.3p as an ordinary differential equation on %. With the nonlinearity defined as 
above it reads 



f £*(T)=L*(r) + N(*(r))forr>-to g r ( , , 

\ $(-logT) = u 1 > 



for $ : [— logT, oo) —> % and initial data u £ %. We rewrite the above system as an integral 
equation, 

$( r ) = S(t + logT)u+ [ 5(r - r , )N($(r , ))dr' for r>-logT, 

J-logT 

i.e., we are looking for mild solutions of Eq. (|4.ip . In order to remove the dependence of the 
equation on the blow up time T we introduce a new variable \£ : [0, oo) — > % defined by 

$r( r ) := $(r - logT) 

such that the above integral equation is now equivalent to 

9(t) = S{t)u+ I S{t - r / )N(^(r / ))dT / for r > 0. (4.2) 

JO 

We study this equation on a Banach space X defined as 



X := <V e C([0,oo),H) : supe^ r ||^(r)|| < oo 

I T>0 

with norm 

H^H* := supe^ T ||*(r)|| 

T>0 

where 

fj, p := \u p \ - e = min{l, _ 1} _ £j 

cf. Proposition 13.91 where e > is arbitrary but fixed and without loss of generality we assume 
e so small that fj, p > 0. In the following, estimate (|3.9p will be used frequently, hence most of 
the constants will depend on e. However, for notational convenience we will only indicate this 
dependence in the proof of the main result. 
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4.4. Global existence for corrected (small) initial data. We follow the strategy of |6J. First 
we study the following equation, 

*(r) = S(t)(1-P)u- / e T - T 'PN(^(T'))dT' + S(t - r , )N(^(r / ))dr / . (4.3) 

Jo Jo 

This is the original equation modified by a correction term in order to suppress the instability 
coming from the symmetry mode. We use a fixed point argument to show existence of solutions of 
Eq. (|4,3p . In a further step we account for the time translation symmetry of the problem and show 
that the correction can be annihilated by adjusting the blow up time T (which is now encoded in 
the initial data) such that we end up with a solution of Eq. (|4.2p . For 5 > we define Xg C X by 

X s := G X : H^ll* < 5}. 

Lemma 4.5. For < 5 < 1 there exist constants c\,C2 > such that 

||N(*(r))|| < Cl 5 2 er 2 ^ T 

and 

||N(tf (r)) - N($(r))|| < c 2 5e-^ T \\y(T) - $(r)|| 

for $,$G^ and r > 0. 

Proof. Let ^ G Xg. Then ||*(r)|| < 5e~^ T < 1 for all r > and 8 < 1. Lemma IP1 implies that 
there exists a constant c\ > such that 

||N(*(r))|| < ci||* (r)|| 2 < Cl 5 2 e- 2 ^ T . 

Let $ G Xg. Then there exists a constant C2 > such that 

||N(*(t)) - N(*(t))|| < |([|*(t)[| + ||*(r)[|)[|*(r) - *(r)[| < c 2 ^HI*(t) - *(t)[| 
which implies the second estimate. □ 
We abbreviate the right hand side of Eq. (14. 3D by defining the operator 

K(*,u)(r) := S(t)(1- P)u- / e r_r 'PN(tf (r^dr 7 + / 5(r - r )N(*(r , ))dT'. (4.4) 

Lemma 4.6. For (5 > sufficiently small and fixed u G H, with ||u|| < <5 2 , the operator K maps 
Af,5 into itself and is contracting, in particular 

||K(*,u)-K(tf,u)||* < ~||*-*||* 

for 3>, \I / G Xg. 

Proof. For fixed u) with ^ G Xg and u G % the integrals occuring in the operator K can 
be viewed as Riemann integrals over continuous functions, which exist since [|PN(\I/(t))|| < 1 by 
Lemma 14.51 To see that K( X I', u) G Xg for 5 small enough we decompose the operator according to 

K(tf , u) = PK(f, u) + (1 - P)K(¥, u). 

We apply the results of Proposition 13.91 and Lemma 14.51 Let ||u|| < 5 2 . Then for r > we obtain 



|PK(tf,u)(r) 



/ e T - T 'p-N(ty(T'))dT' + / S(t - r / )PN(^(r / ))dr / 
Jo Jo 

/OO /'OO 
e T - r '||PN(f {T'))\\dr' < Cl 5 2 / e T - r '( 1+2 ^W < 5 2 e' 2 ^ 
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and 



P)K(*,u)(r)|| < \\S(t)(1-P)u\\+ [ T \\S(T-T')(l-P)lSi(^(r'))\\dr' 

Jo 

< e-^ T ||u||+ [ T e-^ T - T ^\\N(^(T'))\\dT' <5 2 e-^ T + 5 2 f e'^+^dr' 
Jo Jo 



< 5 2 e~^ T . 

We infer that there exist constants c±, C2 > such that 

||PK(*,u)(t)|| < Cl 5 2 e'^ T , 

and 

||(l-P)K(*,u)(r)|| <c 2 6 2 e-^. 

_ L txti 

2ci ' 2c 2 



Thus for 5 < min{l, oT - , oT - } we obtain 



||K(*,u)(r)|| < ||PK(*,u)(r)|| + ||(1 - P)K(*,u)(r)|| < S -e~^ + \e~^ < 5e~^. 

Continuity of K(^, u)(r) as a function of r follows essentially from strong continuity of the semi- 
group (cf. Lemma 3.10 in [6]). It is left to show that K is contracting. Let ^, $ G X$. Then 

e^'HPN^r')) - PN(^(r / ))||dr' 

/"CO /"CO 

< 6 e T - T ' il+ ^\\^T')-V(T')\\dT' <6supe^ a \\^a)-V{a)\\ e T - T ' (1+2 ^W 

Jt cr > T Jt 

< 5e- 2 ^ T \\9-9\\ x . 
Similarly, 

||(1 - F)K(*,u)(r) - (1 - P)K(*,u)(r)|| < f \\S{r - r')(l - P)(N($(r')) - N(*(r')))Mr' 

Jo 

< f e-^(— ')||N($(rO) - N{W))\\dr' < 6 f e^l^r') - ^(r')\\dr' 
Jo Jo 

< 5 sup e^ a \\$(a)-V{a)\\ I e-^ T+T 'W' < 5e~^ T ||$ - ^\\ x . 

o-e(0,r) Jo 

This shows that for 5 sufficiently small, 

supe^||PK($,u)(r)-PK(*,u)(r)|| < h\<f> - y\\ x , 
r>0 4 

and 

supe^Kl - P)K($,u)(r) - (1 - P)K(*, u)(r)|| < j||<3? - *||*, 
r>0 4 

which implies the claim. □ 

Theorem 4.7. For u £ B\ C H sufficiently small, there exists a unique solution ^(-;u) E X of 

*(.;u) = K(tf(.;u),u). (4.5) 

Moreover, the map ^Sf :U C B\ — )• X defined by *(u) = \I>(-; u) is continuous and Frechet differen- 
tiate at u = where U denotes a sufficiently small open neighbourhood of zero in rl. 
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Proof. Lemma 14.61 and the fact that X$ is a closed subset yield a unique fixed point of Eq. (j4.5j) in 
X$. That this is indeed the unique solution in the whole space X follows by standard arguments 
(see also the proof of Theorem 14. lip . Note that for u,v£Wwe have ^(u), \l/(v) E X$ and 

||¥(u)-¥(v)||* = ||K(*(u),u)-K(¥(v),v)||* 
< ||K(*(u), u) - K(*(v), u)||* + ||K(*(v), u) - K(*(v), v)|U- 

By Lemma |4.6[ 

||K(*(u),u)-K(*(v),u)||* < ~||¥(u)-*(v)||*. 
Inserting the definition of K yields 

||K(*(v), u)(r) - K(*(v), v)(r)|| = ||5(r)(l - P)(u - v)|| < e"^||u - v|| 
and we conclude that 

||¥(u)-¥(v)||* < ||u-v||, (4.6) 
which implies continuity. We claim that the solution map \I/ is Frechet differentiable at u = 0. We 
define an auxiliary operator £>*(0) : U -)• X by [.D*(0)v](r) := S(t)(1 - P)v for v € U. It is 
obvious that this defines a bounded linear operator from H into X. We show that it is indeed the 
Frechet derivative, i.e., 

lim JL||*( V ) - *(0) - Z?¥(0)v||* = 0. 

v— >0 v 

Recall that N(0) = 0, hence *(0) = is a solution of Eq. fj4.5j) for u = 0. We assume that v G U, 
such that \I/(v) = K(\I/(v), v). Inserting the definition of K we compute 

*(v)(r) -5(r)(l-P)v= / 5(r - tON^vXt'))^' _ / e — 'pN(*(v)(r / ))dT' 

Jo Jo 
=: G(*(v))(t). 

Again we write G(*(v))(r) = P[G(*(v))(r)] + (1 - P)[G(*(v))(r)]. Estimate g^D and calcula- 
tions similar to those in the proof of Lemma 14.61 yield 

/oo roo 
'||PN(*(v)(r'))Mr'< jf '||*(v)(r')|| 2 dr' 

/>oo 

< || v ||2 / e T-r'(l+2^ p ) dr / < || v ||2 e -2 Mp ^ 



Similarly, 

||(l-P)[G(*(v))(r)]|| < r||5(r-r')(l-P)N(*(v)(r'))Mr' 

JO 

< r e -M P (r-r ')|| N (*(v)(^))||d^ < ||v||V 

We infer that || G( 1 I r (v)) ||^t- ^ ||v|| 2 and thus lim v ^ jrar||G ! (*( v ))||# = 0, which implies the claim. 

□ 

4.5. Global existence for arbitrary (small) initial data. The aim of this section is to use 
the existence result of Theorem 14.71 to obtain a solution of the original wave equation for arbitrary 
initial data (close to 4> T )- Up to now we implicitly assumed the blow up time T to be fixed. 
However, arbitrary perturbations of the initial data will change the blow up time and we account 
for this fact by allowing T to vary. Recall that the initial data we want to prescribe are of the form 

*(o)(p)=I ^ h (47) 

(Tpf'(Tp) + f(Tp))-^- 1 
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see Eq. (|2.3|) . We separate the dependence on T and the free data (/, g) by introducing 



i 



v(p) := I P9{P) P^' 1 ^ , K(p) := ^ (4-8) 

1 pf(p) + f( P ) - K<r / v 1 

which are the initial data relative to the fundamental self-similar solution for T = 1. We rewrite 
the right hand side of (|4.7p and define 

U(v,T)(p) := T^[v(Tp) + «(Tp)] - K {p). 

The data have to be prescribed on the interval [0, T] and we are confronted with the problem that 
we do not know T in advance. As in [6] the argument will be perturbative around T = 1 and 

rl 3^ 

U :=L 2 (0, §) x L 2 (0, |) 



therefore it suffices to restrict T to the interval I = (A, |). In the following we set 



Lemma 4.8. The function U:?ix/->^ is continuous and U(0, 1) = 0. Furthermore U(0, •) : 
I 7i is Frechet differ entiable and 



i 



[D T \J(0,T)\ T=1 X](p) = ^ K »-i g (p), 
where A G R and g denotes the symmetry mode (cf. Eq. ft3.4\ ))- 

Proof. The proof of continuity is similar to the proof of Lemma 3.14 in [6]. We define J : L 2 (0, |) x 

2' 



I -> L 2 (0, 1) by J(v,T)(p) := u(7». For fixed T the map J(-,T) : L 2 (0, |) -> L 2 (0, 1) is Lipschitz- 



continuous since 

\\J(v, T) - J(v, T)|| 2 2(0jl) = jf ^(rp) - v(Tp)\ 2 dp = ±£ \v(p) - v(p)\ 2 dp < 2\\v - £|| 2 2(0if) 

and the continuity is uniform with respect to T. It is therefore sufficient to show that for fixed 
v G L 2 (0, |) the function J(v, ■):/—> £ 2 (0, 1) is continuous. This can be seen by noting that for 
all v,v G L 2 (0, |) and T,T G / 

\\J(v,T) - J(v,f)\\ L 2 {0)1) < \\J(v,T) - J(v,T)\\ L2m + \\J(v,T) - J(v,f)\\ L 2 m 

+ \\J(v,f)-J( V ,f)\\ L 2 m 

i$ \\ v - ^11^(0,1) + ll J (^ T ) - J (v,f)\\ L 2 (0A) . 
Thus, for any given e > we can find a v G C[0, I] such that 



||J(« 1 T)-J(«,r)|| ra(0il) <-+c / |5(Tp)-£(Tp)| 2 dp (4. 



o 



2 

for some constant c > since C[0, |] is dense in L 2 (0, |). By the continuity of v, the integral 
vanishes in the limit T — > T. The above results imply continuity of J : L 2 (0, |) x / — > L 2 (0, 1) and 
thus, 

U(v T)=( T ^( J K T ) + J(ki,T)) -kA 
\ T^(J(v 2 ,T) + J(k 2 ,T))-k 2 J 

is continuous for v = (vi, V2) T G % and re = («i 3 K2) t as defined in Eq. (|4.8j> . 

To show differentiability we set v = and consider U(0, ■) : J — )■ H, which is given by 



U(0,T)(p) = T~k{T P ) - k( P ) = k?- 1 ' p- 1 



| 2 

Tp-i - 1 
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The map is obviously differentiable for all T € I. Recalling the definition of the symmetry mode 
in Lemma 13.51 we obtain 



[L> T U(0,T)| T=1 A](p) 



2A 
p — 1 



.p-i 




2A 

/ 

p — 1 



.p-i 



s(p) 



for A G 



□ 



With these technical results at hand we now turn to the original problem. In the previous section 
we showed existence of solutions for the modified integral equation (|4,3p with initial data u G U, 
where IA denotes a sufficiently small neighbourhood of G H. We rewrite the initial data in terms 
of T and v as defined in Eq. (|4.8p . Inserting in the definition yields U(0, 1) = 0. By continuity 
U(v,T) G U provided that (v,T) G V x / where V and / are sufficiently small neighbourhoods of 
G H and 1 G /, respectively. By Theorem 14. 71 there exists a solution U(v, T) \— > \l/(U(v, T)) G X. 
Recall that Eq. (|4.3|) is Eq. (I4.2p modified by an expontential factor times the function F : V x / — > 
(g) defined by 

F(v,T) :=P(U(V,T)+ f°° e -'N(*(U(v, T))(r'))dr' 



Evaluation yields F(0, 1) = 0, i.e., for v = and T = 1 the correction vanishes and *(U(0, 1)) = 
is also a solution of Eq. (|4.2p . In the following we show that for every small v there exists a T close 
to one, such that this still holds true. We need the next lemma as a prerequisite. 



Lemma 4.9. F: VxIcHxI 
differentiable at T = 1 and 



(g) is continuous. Moreover F(0, 



I 



(g) is Frechet 



D T F(0,T)| T=1 A = ^j 



2A 



/or A G 



Proof. To rewrite F in a more abstract way we introduce the integral operator B : X — > H, \l/ i— > 
Jo°° e ~ T x ^( r/ )^ r/ > which is linear and bounded since 



[•CO 

|B*|| < / e~ T ' ||*(r') Hdr 7 < sup ||*(t')|| < ||^| 

JO r'>0 



A'- 



We define N : X -)■ <Y by N(*)(r) := N(*(r)). We claim that N is Frechet differentiable at £ X 
and the Frechet derivative at zero is given by DN(0)*I' = for ^ G X. This follow from N(0) = 
and 

||N(\I/) \\% ../'(•"! TVT/>T</_>M «.'* ../'(."i >Ti ( -\ i - IUT<i - 



supe^ r ||N(^(r))|| < supe^ T ||*(r)|| 2 < ||*||^ for * G X s . 

T>0 T>0 



Thus 



|N(*)II 



a- 



< 11*1 



x 



X 



with a constant independent of *, which implies the claim. Now 

F(v, T) = P [u(v, T) + BN(*(U(v, T))) 

By Lemma 14.81 and the continuity of N and respectively, we see that F is continuous. To show 
differentiability we set v = and obtain 

F(0, T) = P [U(0, T) + BN(*(U(0, T))) 
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The right hand side is differentiable at T = 1 by Theorem 14.71 Lemma 14.81 and the above consider- 
ations. We conclude that 



L> T F(0,T)| T=1 A = PD T V(0,T)\ T=1 X + PBDN{0)D^(0)D T V(0,T)\ T=1 X 
= PZ? T U(0,T)| T=1 A = ^ r ^g. 

□ 

Lemma 4.10. Let V C H be a sufficiently small neighbourhood ofO. For every v G V there exists 
aT G I C (|, |), such that F(v,T) = 0. 

Proof. The range of F is contained in (g), which is a one dimensional vector space. Thus, there exists 
an isomorphism i : (g) —> R such that i(cg) = c for c£l. We set / := i o F, where / : V x / — > R 
is continuous and F(0, 1) = implies /(0, 1) = 0. Lemma 14.91 shows that /(0, •) : I — > R is 
differentiable at T = 1 and Dr/(0, T)| T=1 7^ 0. Consequently, there exist values T\,T% G I such 
that /(0, Ti) > and /(0, T2) < 0. Continuity of / with respect to the first variable implies that 
there exists an open neighbourhood V C V such that /(v, T\) > and /(v, T2) < for v G V. For 
v G V consider /(v, •) : / — > R. By continuity of /(v, T) with respect to T and the intermediate 
value theorem we conclude that there exists a T* G (7\,T2) such that /(v, T*) = 0. □ 

This yields the next result. 

Theorem 4.11. Let v G H be sufficiently small. Then there exists a T close to 1 such that 

*( T ) = 5(r)U(v,T) + / S(t- r , )N(*(r / ))(ir / , r > (4.10) 

has a continuous solution ^ : [0, 00) — > T~L satisfying 

||*(r)|| < e5e^ r 

/or all t > and some 5 G (0, 1). Moreover, this solution is unique in C([0, 00), %). 

Proof. The existence of a solution \I/ G X$ follows from the above considerations. Let <3? G 
C([0, 00),%) be another solution satisfying the same equation. We assume that ^ 7^ By 
continuity, there exists an e G (0, and a to > such that 

e < ||*(r ) - $(t )|| 

and 

||*(r)-$(r)|| <2e, tG[0,t ], 
which yields ||$(r)|| < 1. For r G [0,To] we obtain 

||*(t)-*(t)|| < c [ T e T - T '\\N(^(T'))-N(^(T'))\\dT' 
Jo 

< C(r )(e T -l) sup ||tf(r') - $(r')|| 

r'e[0,r] 



by applying Lemma 14.41 We infer that there exists a ti G (0, To] such that 

sup ||*(t) - $(t)|| < - sup ||^(t) - $(r)|| 
re[o,n] 2 re [o,n] 

which implies ^(t) = ^(t) for all r G [0,ti]. Iterating this argument yields ^(t) = ^(t) for 
r G [0, t ], which contradicts ||*(r ) - $(t )|| > e. □ 
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Proposition 4.12. (Global existence for arbitrary, small initial data) Let e > be small enough 
such that p p = \oj p \ — e > 0. Let v 6 L 2 (0, |) x L 2 (0, |) 6e sufficiently small. Then there exists a 
T close to 1 such that 

$( T ) = 5(r + logT)U(v,T) + / 5(r - r / )N($(r , ))rfr / , r>-logT (4.11) 

./-logT 

/ios a continuous solution : [— log T, oo) — > T~L satisfying 

\Mt)\\ <C e e~^ T 

for all t > — logT and a constant C £ > depending on e. Moreover, this solution is unique 
in C([— logT, oo), ■%). Thus, $ is i/ie unique global mild solution of Eq. 1\ ) with initial data 
<&(-logT) = U(v,T). 

4.6. Proof of Theorem [Ql 

Proof. We translate the result of Proposition 14. 121 back to the original coordinates (t, r). Let (/, #) 
satisfy the assumption of Theorem 11.21 For the fundamental self-similar solution with T = 1 we 
have 



We define 



^ 1 (0,r) = « *" 1 . ^ 1 (0,r) = ^ I K p - 1 



:= - -^r^- 1 , v 2 (p) := /(p) + pf'(p) - kE" 1 



v 

such that v = (vi,v 2 ) T G L 2 (0, |) x L 2 (0, |) and 



2>' -^ V> 2, 

H« = ||(/^)-(^ 1 (o, 0,^(0, OJIUd) 

We may assume v small enough to satisfy the assumptions of Proposition 14.121 and we infer that 
there exists a unique global mild solution $ € C([— logT, oo),%) of Eq. (|4.ip for T close to 1 with 
initial data 3>(— logT) = U(v, T) and 

||*(r)|| < C £ e~^- £ ^ T 

for all r > — log T. By definition 

d>(r)(p) = (0 1 (r,p),0 2 (r, /O )f 
is a solution of Eq. (pOI) and Eq. ([23]) yields 

V>(i, r) = ^ T (t, r) + (T - t)~^r~ 1 / 2 (- log(T - i), ^)dr' 



and 

r) = ${t, r) + (T - t)-^r r - log(T - t), ^). 
For ip = ip — we obtain 



||fo(t,.), ^(*,0ll!(T-t) = (T-t)-^r( / T V2(-log(T-t),^)| 2 dr+ T log(T - i), t^)| 2 ^ 



= (T-t)—^ \M-log(T-t),p)\'dp + ^ \M-log(T-t),p)\'dp 

= (T -t)^\\<S>(-\og(T -t))\\ 2 < C e 2 (T-i)^T +2(l ^ | - E) . 
Thus, 

||ty(t, V*(*> •)) " ty T (*> 0. ^ O)IU(r-t) < C £ (T - t )^) + ^l- £ . 

□ 
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